Collective effects in spin-crossover chains with exchange interaction 
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The collective properties of spin crossover chains are studied. Spin crossover compounds contain 
ions with a low-spin ground state and low lying high-spin excited states and are of interest for 
molecular memory applications. Some of them naturally form one-dimensional chains. Elastic 
interaction and Ising exchange interaction are taken into account. The transfer-matrix approach is 
used to calculate the partition function, the fraction of ions in the high-spin state, the magnetization, 
susceptibility etc. exactly. The high-spin/low-spin degree of freedom leads to collective effects not 
present in simple spin chains. The ground state phase diagram is mapped out and compared to the 
case with Heisenberg exchange interaction. The various phases give rise to characteristic behavior 
at nonzero temperatures, including sharp crossovers between low- and high-temperature regimes. 
A Curie- Weiss law for the susceptibility is derived and the paramagnetic Curie temperature is 
calculated. Possible experiments to determine the exchange coupling are discussed. 

PACS numbers: 75.10.Hk, 75.20.Ck, 75.50.Xx 
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I. INTRODUCTION 

Motivated in part by the search for molecular memory 
deviceSfiiSi^ spin-crossover compounds (SCCs) have been 
investigated quite intensively in recent years. These com- 
pounds are characterized by magnetic ions that can be 
either in a low-spin (LS) or high-spin (HS) state, where 
the HS state is slightly higher in energy and can be ther- 
mally populated4*^*SiL& An ion in the HS (LS) state has 
the spin quantum number 5'hs (5'ls), where S'hs > 'S'ls. 
Most SCCs consist of metal-organic complexes involv- 
ing transition- metal ions, in many cases Fe^"*". Spin 
crossover involving charge transfer between two differ- 
ent ions is observed in Prussian Blue analogues i2ii2iii 
SCCs are promising for molecular memory applications 
due to long lifetimes of the HS state. Applications re- 
quire to deposit SCCs on substrates as thin films or one- 
dimensional (ID) chains. Many SCCs consist of weakly 
interacting chains or plains even in the bulki^iiSiiii^ A 
typical quasi-lD material is Fe^"*" with 4-R-l,2,4-triazole 
ligandsi^^ 

SCCs are also interesting from the point of view of 
basic physics. The additional LS/HS degree of freedom 
leads to collective effects not present in pure spin models, 
which are most pronounced for a diamagnetic LS state, 
•S'ls = 0, as in the case of Fe^^ ions. Then, in the LS 
state the spin is switched off. This reminds one of diluted 
spin modelsjiSiiSiil but with fluctuating dilution. 

Most of the theoretical literature on SCCs omits the 
exchange interaction. The only interaction in this case 
is of elastic origin and determines whether neighboring 
ions prefer to be in the same (LS or HS) state or in dif- 
ferent onesii^ii^ However, typical values for the exchange 
interaction J in transition-metal complex salts are of the 
order of |J|/fcB = 10 to 20Ki^W^ The interaction is 
due to superexchange and is typically antiferromagnetic. 
In ferrimagnetic Prussian Blue analogues the exchange 
interaction is typically similariSiifl but values of the order 



of room temperature have also been reported^ 

If magnetic anisotropies are small, the exchange in- 
teraction is of Heisenberg type. For strong easy-axis 
anisotropy we can instead consider an Ising exchange in- 
teraction. For isolated ID chains deposited on a sur- 
face strong easy-axis anisotropy is expected due to shape 
anisotropy. Another important source of anisotropy is 
the ligand field acting on the magnetic ions. Nishino 
et a;,23i24,25 g^jjj Boukheddaden et alM. study a three- 
dimensional model with nearest-neighbor Ising exchange, 
S'ls = 0, and S'hs = 1/2 (the high spins have two 
possible orientations) employing mean-field theory and 
Monte Carlo simulations. On the other hand, Timm and 
SchoUwock^^ consider a ID model with nearest-neighbor 
Heisenberg exchange, Sls = 0, and Shs = 2. The zero- 
temperature phase diagram is investigated using the den- 
sity matrix renormalization group:-- 

In this paper, we study the behavior of a ID spin- 
crossover chain with and without Ising exchange interac- 
tion at all temperatures. We will argue that its solution 
also provides a good approximation for the Heisenberg 
case for small exchange interactions. In Sec.^we present 
the model and discuss the transfer-matrix calculations. 
We mostly restrict ourselves to the case Sls = and 
Shs = 2, which is appropriate for Fe^^ compounds. The 
results are presented in Sees. IIIll and II VI for the ground- 
state and finite-temperature properties, respectively. 



II. MODEL AND METHOD 

This section outlines the theory for general Shs: unless 
stated otherwise. We start from the Hamiltonian 



H 
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where = with rri; = for (t.; = 1 (LS state) and 
rrii = — Shs,---iShs for = ^1 (HS state). Here, 
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V describes the elastic interaction between neighboring 
spins>i^ For V > {V < 0) homogeneous (alternating) 
arrangements of LS and HS are favoredii^ 2A is the en- 
ergy difference between HS and LS and h is the magnetic 
field with a factor gjiB absorbed. J > (J < 0) is a fer- 
romagnetic (antiferromagnetic) exchange coupling. For 
J = this model is equivalent to the one of Ref. 29 and 
for J = /i = it is the model introduced by Wajnflasz 
and Pick^ and by Doniachi^ 

We apply the transfer-matrix approach. To focus on 
the essential physics we assume that the vibrational fre- 
quencies of a complex do not change between the LS and 
HS state. Then the degeneracy of the HS state is only 
due to its spin. The main effects of different vibrational 
frequencies would be to renormalize the energy splitting 
2A and to change the effective degeneracies of the LS and 
HS statesiSiSS These effects do not change the qualitative 
results and can easily be reintroduced. 



The total partition function is Z = ^ 



where /3 
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^ ^ (2) 



The partition function then reads 

^ — ^ ^ g-f^cri7Tii,cr2Tn2 g-^cr2TrL2, 0-3771 

{(Ti ,mi } 

(3) 

This can be written as Z ~ Tr where the sym- 
metric matrix M has the components Mcm,a'm' — 
eicp{Ka-m,.tT'm')- For iV — !■ 00 the partition function be- 
comes the maximum eigenvalue of M to the power N and 
the Gibbs' free energy per site is 

g = -T lim — InTrAf^. (4) 



For J ~ the eigenvalue equation can be solved in closed 
form. For 5ls = we obtain 
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= -i In ( e'^^ cosh PA + J e^P^ sinh^ /3A 



/3 



-2I3V 



In .gns 

2/3 



(5) 



with A = A - ln5Hs/2/3 andSa 

g/9/i(SHS + l/2) _ g-/3h(SHS + l/2) 



5HS 



^0hm 



m=-S'HS 



-ph/2 



(6) 

For h = Q one has gns = 2S'hs + 1- 

For J ^ 0, ^Ls = 0, ^Hs = 2 the calculation of the 
partition function has been reduced to the eigenvalue 
problem of the 6x6 matrix M. An important quan- 
tity describing SCCs is the fraction 7 of ions in the HS 
state. It is given by 7 = (1 — (o'))/2 so that, with Eq. ||SJ), 
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Tr Af ^ ■ 
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FIG. 1: (Color online) Ground-state phase diagram of the 
spin-crossover chain without exchange interaction and with 
Sls = 0, in terms of elastic interaction V and Zeeman energy 
h in units of A, where 2 A is the LS/HS energy splitting. The 
solid lines denote discontinuous transitions. 



If |n), n = 1, . . . , 6 are the orthonormalized eigenvectors 
of M with eigenvalues mi > m2 > • • • > me, then 

1 T (1|9A//9A|1) 
2 
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2 mi 
Similarly, the magnetization is 

dg ^ {l\dM/dh\l) 



(8) 



mi 



(9) 



and the probability of any two neighbors being in the 
same (LS or HS) state is 



1-dg/dV 



T {l\dM/dV\l) 
2" 



mi 



(10) 



This quantity describes nearest-neighbor correlations. 
For the special case J = we find the HS fraction 
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s'^^ sinh f3A 



;2/3V sinh^ /3A + e 



-2I3V 



(11) 



In this case the magnetization is determined by 7 through 
(m) = 75'hs-Bshs (/^^'S'hs)j where Bs(x) is the Brillouin 
function. This leads to a simple result for the suscepti- 
bility X = d{m)/dh for vanishing magnetic field: 



X = 7- 



>5'hs('5'] 



HS 



1) 



3T 



(12) 



Thus the susceptibility is just the paramagnetic expres- 
sion weighted by the concentration of high spins. 



III. GROUND STATES 

We now derive the ground-state phase diagram to lay 
the ground for the discussion of finite-temperature prop- 
erties. For the case without exchange interaction Fig. 
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FIG. 2: (Color online) Ground-state phase diagram of the 
spin-crossover chain with Ising exchange coupling J. In this 
and the following graphs Sls = and S'hs = 2 are assumed. 
The planes denote discontinuous transitions. The heavy solid 
lines show their intersections with the J — plane, cf. Fig. 
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FIG. 3: (Color online) Solid lines: Intersection of the T = 
phase diagram, Fig. |21 with the V = plane (vanish- 
ing elastic interaction). Dashed curves: The phase bound- 
aries for Heisenberg exchange interaction!— Long-dashed 
(short-dashed) curves correspond to discontinuous (contin- 
uous) transitions. The alternating phase is reduced to two 
finite lobes by quantum effects. Note that the phase bound- 
aries for the Ising- and Heisenberg- type models are identical 
for positive and small negative J. 



which does not exist for J = 0. It becomes unstable at 



shows the phase diagram in terms of the two dimcn- 
sionless ratios V/IS. and ft./ A. From Eq. Q we find 
7 = 0, (m) = 0, and Wcq = 1 for |/i|/A < 2/5'hs and 
F/A > -1/2 + |ft|5Hs/4A. This is the LS state. For 
|/i|/A > 2/5'hs and F/A > 1/2 - |/i|S'hs/4A we find 
7=1, (m) = S'nssgn/i, and Weq = 1- This is the fully 
polarized HS state. In all other cases we find 7 = 1/2, 
(m) = S'HsSgnft/2, and = 0, corresponding to an 
alternating state of LS and HS ions. 

In the general case with Ising exchange interaction we 
have three dimensionless parameters, 1^/A, J/A, and 
h/ IS.. The ground state is found by determining the 
largest component(s) of the matrix M, since for T ^ 
these components become exponentially larger than the 
others. For sufficiently large J/A > we obtain a ferro- 
magnetically aligned HS phase for any V, h. For ^hs = 2 
and h> this phase becomes unstable at 
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(13) 

For large — J/A > we find a HS phase with Neel order, 
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(14) 

Finally, for small |J/A| the LS and alternating phases 
are separated by a boundary at 
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A 
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2A 



for 



1 J 1 

2<A<2 
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(15) 



The resulting T = phase diagram is shown in Fig. |21 
Figure |3| shows the V — Q cut of the phase diagram and 
compares the results to the case of Heisenberg exchange 
interaction from Ref. In the Heisenberg-type model 
several more complex phases are found, which are absent 
in the Ising case. These consist of arrangements of 2, 3, 
5, 7, or 9 consecutive ions in the HS state followed by 
a single ion in the LS state. This pattern is periodically 
repeated. These phases are stabilized in the Heisenberg 
case by the energy gain for antiferromagnetically coupled 
spins due to quantum effects. This energy gain also sta- 
bilizes the antiferromagnetically aligned HS state, which 
displaces the other phases, as seen in Fig. |3| 

On the other hand, for all positive and small negative 
J/A the phase boundaries coincide. This is because the 
LS, alternating, and ferromagnetic HS ground states of 
the Ising-type model remain eigenstates, with the same 
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energy, of the Heisenberg-type model. (The LS and alter- 
nating states do not contain nearest-neighbor HS pairs so 
that the exchange interaction does not enter. The ferro- 
magnetically fully aligned state is an eigenstate for any — 
even antiferromagnetic — Heisenberg-type model.) In ad- 
dition, all ground states that appear at larger |J| have 
nonzero energy gaps to the ground state at small |J|. 
Thus for small |J|, which are expected for most SCCs, 
no new phases appear for either the Ising or the Heisen- 
berg case and both have the same ground states. 



IV. FINITE-TEMPERATURE BEHAVIOR 

We now turn to the finite-temperature properties. The 
partition function is analytic for T > so that the phase 
transitions are replaced by crossovers. These can become 
very sharp, however. 

While the assumption of Ising exchange is appropriate 
in the presence of strong easy-axis anisotropy, it also pro- 
vides a reasonable approximation for the Heisenberg case 
for small J. In this case the ground state is found ex- 
actly and the energies of low-lying excited states above 
the LS and alternating ground states are also identical 
for both cases so that the behavior at low temperatures 
will be very similar. Deviations will appear when a sig- 
nificant fraction of nearest neighbors are both in the HS 
state since then the exchange interaction becomes rele- 
vant. We will come back to this at the end of this section. 
On the other hand, in the ferromagnetic HS state the low- 
energy excitations are different but gapped in both cases. 
Thus only qualitatively similar behavior is expected. 

Note that for solving the Heisenberg-type model, map- 
ping onto and solving an Ising-type model is superior to 
the mean-field approximation, which neglects all fluctua- 
tions. In the Ising approximation we treat fluctuations of 
CTi and longitudinal fluctuations of Si exactly and neglect 
only transverse fluctuations of S^. 
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FIG. 4: (Color online) Solid curves: HS fraction 7 as a func- 
tion of temperature for various values of the elastic interaction 
V and J — fi — 0. Dashed curves: 7 for h/ A — 0.8 and J = 0. 
The ground state is in the LS phase. 
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FIG. 5: (Color online) Solid curves: Magnetization (m) /Sns 
as a function of temperature for various values of the elas- 
tic interaction V, Zeeman energy h/A = 0.8, and exchange 
J — 0. Dashed curves of same color: Probability lUeq of neigh- 
boring ions being in the same (LS or HS) state for the same 
parameters. The ground state is in the LS phase. 



A. Low-spin phase 

Typical SCCs are in the LS phase for T — > 0. Figure 
0] shows the HS fraction as a function of temperature for 
J = and various elastic interaction strengths V. For 
V = we observe the well-known smooth crossover to 
the HS state, which comes from its higher degeneracy. 
For T 00 the HS fraction is only determined by the 
degeneracies, 7 (2S'hs + 1)/(25'hs + 2) = 5/6, for any 
V. For V > the interaction favors a homogeneous state 
and the crossover becomes sharper. In mean-field ttieory 
we would eventually reach a critical point V = Vc and 
find a discontinuous transition for V > Vc^ However, the 
exact solution for ID shows that 7 develops a discontinu- 
ity only for V —f 00. For larger V/ A the curves overshoot 
and exhibit a maximum. This is due to the energy gain 
from the homogeneous HS state outweighing the entropy 
gain from populating all states equally. For V < the 



interaction favors a mixed state with 7 « 1/2 and thus 
broadens the crossover. 

All curves for fixed h in Fig. 0] cross at 7 = 1/2. The 
equation 7(Tly) = 1/2 is equivalent to (a) = and, using 
Eq. ||ni), to A - (T/2) Ingns = 0. This leads to the im- 
plicit equation T^/A = 2/lTigiis{h,T^), which is indeed 
independent of the interaction V so that all curves cross 
at T^. For /i = we obtain T^/A = 2/ ln(25Hs + !)■ 

The magnetization (to) and probability Wcq for equal 
neighbors are shown in Fig. Ofor h/A — 0.8. For V = 
we see a broad crossover in (to) from the LS phase to a 
maximum and then a decay ~ 1/T at high temperatures. 
The Curie- Weiss law is discussed further in Sec. lIV CI For 
V > {V <0) the crossover becomes sharper (broader). 
Similar behavior is found for J ^ 0. The nonmonotonic 
temperature dependence can be interpreted as the ID 
analogue of the reentrant transition to ferromagnetic or- 
der found in three dimensionSi^ 
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FIG. 6: (Color online) HS fraction 7 as a function of tem- 
perature for various values of the exchange interaction J and 
V = /i = 0. The ground states are always in the LS phase. 
Inset: HS fraction 7 as a function of exchange interaction 
J/ A for various values of the Zeeman energy h for T/A = 1 
and V = ^. The T = transition between LS and antifer- 
romagnetic HS phases it at J/ A = —1/2, the positions of 
the transitions to the ferromagnetic HS phase are marked by 
triangles of the same color as the corresponding curves. 

The probability Woq starts from Wcq = 1 at T = 
and approaches a universal value of 13/18 for T — s- 00. 
Weq shows a broad minimum for small V ^ which is best 
understood in the case without interactions. Then = 
7^ + (1 — 7)^ so that 7 = 1/2 impKes Wcq = 1/2, which 
is smaller than the limits for low and high T . 

Next, we consider the effect of the exchange interac- 
tion, see Fig. El For vanishing magnetic field the effect 
of moderate | J/A| is rather weak due to the scarcity of 
nearest-neighbor HS pairs. For J 7^ the crossover tem- 
perature with 7(T-y) = 1/2 depends on both J and V . 7 
increases symmetrically for ferromagnetic and antiferro- 
magnetic J for all V ^ due to the invariance of the Hamil- 
tonian under — > (— l)*mi and J — > — J for /i = 0. The 
T = transitions to the ferromagnetic and antiferromag- 
netic HS phases take place at J/A = ±1/2, respectively. 
For /i 7^ the transition to the ferromagnetic phase shifts 
to lower values of J, whereas the antiferromagnetic tran- 
sition remains fixed, see Fig. |21 so that the minimum of 
7(J) shifts to lower J, as seen in the inset of Fig.|H| 

B. Other phases 

The alternating LS/HS phase is stabilized by a strongly 
negative elastic interaction V This can be expected for 
SCCs grown on a substrate, where the transition of an 
ion to the HS state increases its radius, which should 
favor a smaller ionic radius (LS state) of its neighbors. 

Figure |7| shows the HS fraction and the probability 
of neighbors in the same state for J = /i = 0. 7 
first decreases with temperature and shows a universal 
crossing at 7 = 1/2 and T = T-y as in the LS phase 



FIG. 7: (Color online) HS fraction 7 (solid curves) and prob- 
ability uicq of neighbors in the same state (dashed curves of 
same color) as functions of temperature for various values of 
the elastic interaction V and J = h — 0. The ground state is 
in the alternating phase. 




FIG. 8: (Color online) HS fraction 7 (solid curves) and mag- 
netization (m) (dashed curves of same color) as functions of 
temperature for various values of the elastic interaction V, 
J = 0, and h/A = 1.5. The state at T = is in the HS phase. 



[T^/A = 2/ln(2S'HS + 1) for h = 0]. The HS fraction 
first decreases since the alternating order is partially de- 
stroyed {wcq increases) and the energy of spins having 
one LS and one HS neighbor is only determined by A, 
which favors the LS state. Since the effect of J 7^ is 
weak as long as the ground state is not changed, we do 
not show the J dependence here. At zero field 7 is an 
even function of J, as for the LS phase. 

While the exchange interaction has only a weak effect 
in the LS and alternating phases, the situation changes 
for the HS phases. The Jerromagnetically aligned phase 
appears at (very strong) magnetic fields and for strong 
ferromagnetic exchange interaction — both might not be 
realized in SCCs. Typical results for the field-induced fer- 
romagnetic HS state are shown in Fig.|Sl 7 goes through 
a minimum since the ground state of a single ion would be 
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FIG. 9: (Color online) (a) HS fraction 7 (solid curves) and 
magnetization (m) (dashed curves of same color) as functions 
of temperature for various values of the exchange interaction 
J, V = 0, and smaU /i/A = 10"^. (b) The same for fixed 
J/A = 1 and various Zeeman energies. The state at T = is 
in the ferromagnetic HS phase. 



the lowest state of the HS quintet, while the first exited 
state is the LS singlet. Therefore for increasing temper- 
ature the LS state contributes first. The magnetization 
for J = is determined by 7 and a Brillouin function 
containing the external field, as noted in Sec. [H] 

For the ferromagnetic phase induced by strong ex- 
change at small field, Fig. ^ shows typical results for 
the HS fraction and magnetization. 7 and (m) now de- 
pend strongly on J. The HS fraction behaves similarly 
to the case of J = and large field but now the simple 
relation between 7 and (m) no longer holds — it would 
predict a tiny magnetization at h/A = 10~^. Instead, 
the magnetization shows a rather sharp crossover from 
nearly full polarization to low magnetization. Figure|5fb) 
shows the dependance of 7 and (m) on magnetic field. At 
small h/ A ^ 10~^ the HS fraction hardly depends on h, 
but the magnetization shows a strong field dependence. 
There is a strong tendency towards ferromagnetic order 
even at very small h, whereas the magnetization vanishes 
for h = 0. We return to this crossover in Sec. IIVDI 

Finally, the antiferromagnetically aligned HS phase 
only exists for strong antiferromagnetic exchange J/A < 
— 1/2. Figure [TUI shows the HS fraction and the magneti- 
zation for fixed J/A = — 1 in a magnetic field. 7 behaves 
very similarly to the ferromagnetic phase. Fig. |S1 The 
field induces a nonzero magnetization only for T > 
because of the energy gap for spin flips. For Heisenberg 
exchange the temperature dependence would be quite dif- 
ferent: The ground state in that case still has 7=1 but 
no long-range order due to transverse spin fluctuations. 
However, it is also gapped since Sus is an integer n^Si 
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FIG. 10: (Color online) High-spin fraction 7 (solid curves) 
and magnetization (m) (dashed curves of same color, scaled 
by 10) as functions of temperature for various values of the 
elastic interaction V, h/A — 0.1, and J/A = —1. The state 
at T = is in the antiferromagnetic HS phase. 



C. High-temperature limit 

In the present subsection we consider the high-tempe- 
rature expansion of thermodynamic quantities. We dis- 
cuss the Curie- Weiss law for the susceptibility and pos- 
sible experiments to determine the exchange coupling. 

In Sec. |n] the partition function z for a single 
site has been obtained as the largest eigenvalue of 
the matrix M. We write M = M^c + AM with 
(Moo)crm,cr'm' = Hmr^oo Mcr,n,a'7n' = 1 and Calculate 
z perturbatively for small AM. Moo has an eigenstate 
|l)o = (1, 1, 1, 1, 1, 1)/V6 with eigenvalue 2S'hs + 2 = 6 
and all other eigenvalues vanish. Using standard pertur- 
bation theory for z and the corresponding eigenstate |1) 
and expanding in 1 /T we obtain 



z = 6 



8V - 12A 



3T 

1211/2 - leSFA + 81A2 + 225J2 + ISS/i^ 
27r2 

{U8V^ - 846y2A + nWA^ + 2925yj2 
+ lOSOy/i^ - 162A^ - 2700A J2 - UlbAh^ 
+ 4050 J/i^)/ (243T^) +C'(l/r4). (16) 



The Gibbs' free energy per site is then 
the susceptibility is 



-Tlnz and 
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d'^g _ 5 5(41/- 3A + 30J) 



+ 0(1/T^). (17) 



dh^ 3T 27 r2 

Expanding the Curie- Weiss law x — C /{T — 0) we flnd 
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with the paramagnetic Curie temperature 
41/-3A + 30J 
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FIG. 11: (Color online) Inverse susceptibility 1/x as a func- 
tion of temperature for vanishing magnetic field, h = 0. (a) 
For J — and various values of V/ A. (b) For V — and var- 
ious values of J/ A. The solid curves show the exact results, 
while the dashed curves of the same color give the Curie- Weiss 
law X = 5/(3(T - e)) with 9 given by Eq. llTUll . 



Figure [TTl shows the approach of the exact result for 1/x 
to the Curie- Weiss form. 

We see that the high-temperature behavior of the sus- 
ceptibiUty not only depends on the exchange interaction 
J but also on the elastic interaction and the LS/HS 
energy sphtting. Interestingly, the paramagnetic Curie 
temperature Q is nonzero even without any interaction, 
in which case 6 = —A/3. Thus there is a deviation from 
the Curie law for the susceptibihty of noninteracting spin- 
crossover ions. The reason is that the ground state is dia- 
magnetic and thermal activation with T ~ A is required 
to make an ion paramagnetic. Obviously, O < does not 
imply an antiferromagnetic ground state. Similarly, the 
ground state is not always ferromagnetic if O > — the 
tendency towards ferromagnetism can be preempted by 
the crossover to the LS state. This is seen for J/A = 0.2 
and 0.4 in Fig.HIIb). 

Since it is difficult to infer the value or even the sign of 
J from the paramagnetic Curie temperature, we propose 
another route to extract J. For J = h = we have 
X = 27/T, cf. Eq. injl. This identity does not hold for 
nonzero J so that the deviation from it might be a good 
measure of J. We thus consider the ratio /i = xT/2'^ at 
vanishing magnetic field. From Eq. (|16|l we obtain 
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1 + dg/dA 



6 



5(4y-3A) 
54T 



-o(i/r2) (20) 



and, at /i = 0, 

-Td^g/dK^ 



1 + dg/dA 



l + ^+0{l/T'). (21) 



FIG. 12: (Color online) Crossover temperatures for J — 0. 
The r = plane corresponds to Fig. Q The surfaces indicate 
solutions of 7 = 1/2 and Woq = 1/2, respectively. 



The leading high-temperature behavior of ^ is entirely 
determined by J. This implies that the effect of V and A 
on the high-temperature susceptibility only comes from 
the HS fraction. Numerical results suggest that /i > 1 
(/i < 1) for J > (J < 0) at all temperatures, x can 
be measured by standard methods, while the HS fraction 
can be obtained independently from optical transmission 
experiments^ or x-ray-absorption near-edge structure!^ 
(XANES). The first method makes use of the change in 
electronic structure between the LS and HS states, while 
the second relies on the change of atomic distances. Since 
X and 7 are measurable, /i is a promising quantity for the 
determination of the exchange coupling. 



D. Finite-temperature crossovers 

The finite-temperature crossovers found in the preced- 
ing subsections are now studied in more detail. Regard- 
less of the specific ground state, for high temperatures 
the system approaches a limit with 7 — 5/6, (m) = 0, 
and 



7^ (1 - 7)2 = 13/18. If one starts from the 
LS ground state, 7 has to go from zero to 5/6. Figures 
0] and El suggest that j{Tj) = 1/2 is a good definition for 
a crossover temperature T^. Starting from the alternat- 
ing phase with 7 = 1/2, 7 first drops and then recrosses 
7 — 1/2, at least for small J, see Fig. 13 A measure for 
the crossover from the alternating phase is WcqiT^) — 1/2 
since Weq vanishes in the alternating phase and is close 
to unity in the HS state at high temperatures. 

We first consider the case J = 0. Figure El shows 
the crossover temperatures. From Eq. Hll|) . 7(T^) = 1/2 
is equivalent to A = 0. Tj becomes zero for |/i|/A = 
S/S'hs = 1. a nontrivial solution exists for |/i|/A < 1, 
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FIG. 13: (Color online) Crossover temperatures for V — 0, 
defined by (a) 7 = 1/2, (b) w^q = 1/2, (c) (m)/SHS = 1/2. 
The T — plane corresponds to Fig. El Phase boundaries at 
T = are marked by heavy solid lines. 



i.e., for ground states in the LS and part of the alter- 
nating phase. Note that Tj is independent of V. The 
equation u'oq(T![„) — 1/2 has two solutions (reentrance) if 
y < and the ground state is in the LS phase, of. Fig.|Sl 
whereas it has only one solution if the ground state is in 



the alternating phase and none for > 0. vanishes 
at the phase boundaries of the alternating ground state. 
The LS state is destroyed by temperatures T ^ 2A, the 
energy per ion for going from the LS to the HS state. 
On the other hand, the alternating state is stable up to 
T ^ 2| J|, the energy cost for breaking LS/HS pairs. 

For J 7^ the typical behavior can already be seen for 
V = 0. Figures [rST a.b'l show the crossover temperatures 
defined by 7 = 1/2 and Wcq = 1/2 as functions oi h/A 
and J/ A. This should be compared to the phase diagram 
in Fig. 121 The reentrant behavior in Wgq is seen both in 
Fig. Ua and Fig.lTSIb). 

For the alternating ground state we find that there is 
(at least) one nonzero temperature with 7=1/2 for some 
parameter values, as in the J = case. The parameter 
range where this is the case narrows and vanishes for 
fields above h/A « 2.27, as seen in Fig. llBf a). As noted 
above, we expect significant differences between models 
with Ising and Heisenberg exchange if nearest-neighbor 
pairs in the HS state appear. This is the case if both 
7 > 1/2 and Weq > 1/2, i.e., above the crossovers in 
Figs. I13r a). (b). For the HS ground states there will be 
deviations already at small T since the spin excitations 
are different in the Ising and Heisenberg models. 

The crossover from the ferromagnetic HS ground state 
to the high-temperature HS state with (m) = is char- 
acterized by the condition (m) = S'hs/2. The crossover 
becomes very sharp at low magnetic field, cf. Fig. (51 (In 
the Heisenberg case it might be broadened by transverse 
spin fluctuations.) Figure lTl^ c) shows the crossover tem- 
perature for V ^ 0, which depends logarithmically on 
magnetic field for small h. This result is known from the 
ID Ising model and can be understood from the forma- 
tion energy of droplets of reduced spin polarization. 



V. SUMMARY 

The thermodynamic properties of ID SCCs with nea- 
rest-neighbor elastic and Ising exchange interactions have 
been investigated using the exact transfer-matrix ap- 
proach. This is motivated by SCCs naturally forming 
ID chains, like Fe^^ with 4-R-l,2,4-triazole ligands, and 
by possible artifical ID structures for molecular electron- 
ics. The Ising model is appropriate for strong easy-axis 
magnetic anisotropy, for example shape anistropy in ID 
chains. For simplicity, the ratio of the degeneracies of 
the LS and HS state of an ion is assumed to be only due 
to their spin, corresponding to the assumption of equal 
vibrational frequencies :^*2£ 

The HS fraction, magnetization, and spin susceptibil- 
ity are calculated, as well as the probability of neigh- 
boring ions being in the same (LS or HS) state. These 
quantities allow one to map out the T — phase dia- 
gram. A strong magnetic field or a ferromagnetic ex- 
change interaction can drive the system into a ferromag- 
netic HS ground state, whereas an antiferromagnetic ex- 
change interaction can lead to an antiferromagnetic HS 
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ground state. A negative elastic interaction can drive the 
system into an alternating LS/HS ground state. These 
phases give rise to characteristically different behavior 
at finite temperatures. Nonmonotonic temperature de- 
pendence of HS fraction and magnetization arc ubiqui- 
tous. While there are no phase transitions at T > 0, 
the crossovers between the low-temperature behavior and 
the HS-dominated high-temperature limit can become 
very sharp. Spin crossover chains show ID analogues 
of interesting effects found in the bulk, such as reentrant 
transitionai2& and stripe phases 

The susceptibility follows the Curie- Weiss law at high 
temperatures for all ground-state phases. The paramag- 
netic Curie temperature is found to depend not only on 
the exchange interaction but also on elastic interaction 
and LS/HS energy splitting. It is nonzero even for nonin- 
teracting complexes. The exchange coupling J could be 
determined experimentally from the susceptibility x and 



HS fraction 7 using that the high-temperature behavior 
of the ratio depends only on J. 

The typical energy scale in all results is A, half the 
LS /HS energy splitting for a single ion. It would thus be 
desirable to tune this scale to small values by a suitable 
choice of ligands so that large values of the dimensionless 
parameters V^/A, J/A, and /i/A can be reached. This 
would allow to test the rich temperature-dependent ef- 
fects predicted here. 
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